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ABSTRACT 

We present a quantitative analysis of the time-delay effect on the cosmic 
background radiation (CBR) by static gravitational potential of galaxy clusters. 
This is primarily motivated by growing observational evidence that clusters 
have essentially experienced no-evolution since redshift z ~ 1, indicating that 
the contribution of a time-dependent potential to CBR anisotropy discussed in 
literature could be rather small for the dynamically-relaxed clusters. Using the 
softened isothermal sphere model and the universal density profile for the mass 
distribution of rich clusters, we calculate the CBR anisotropy by the time-delay 
effect and compare it with those generated by the thermal and kinematic S-Z 
effects as well as by the transverse motion of clusters. While it is unlikely that 
the time-delay effect is detectable in the current S-Z measurement because of 
its small amplitude of 10~ 6 -10~ 7 and its achromaticity, it nevertheless leads to 
an uncertainty of ~ 10% in the measurement of the kinematic S-Z effect of 
clusters. Future cosmological application of the peculiar velocity of clusters to 
be measured through the S-Z effect should therefore take this uncertainty into 
account. 

Subject headings: cosmic microwave background — cosmology: theory - 
galaxies: clusters — gravitation 
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1. Introduction 

The microwave sky behind a cluster of galaxies would be primarily distorted through 
the so-called Sunyaev-Zel'dovich (Z-S) effect - the inverse Compton scattering of the 
radiation by electrons in hot intracluster gas, which leads to a temperature decrement of 
typically AT/T ~ 10~ 4 in the Rayleigh- Jeans part of the spectrum for a cluster as rich 
as Coma. The original goal of conducting the S-Z measurement is to estimate the Hubble 
constant H , while the recent attempt has also been made to determine the cluster baryonic 
(gas) mass in combination with X-ray observation (Myers et al. 1997; Mason & Myers 
2000). Yet, the S-Z effect can only probe the distributions of intracluster gas, and provides 
no direct information about the underlying gravitational potential of the clusters, whereas 
the latter plays a much more important role in cosmological study. This motivates us to 
address the following question: Can we simultaneously detect the gravitational imprint on 
the microwave sky by cluster potential during the measurement of the S-Z effect ? 

Essentially, the gravitational potential of a cluster exercises an influence on the cosmic 
background radiation (CBR) through the gravitational lensing effect and the Rees-Sciama 
effect (Rees & Sciama 1968), which consist primarily of three components: The first one 
arises from gravitational lensing, which alters the trajectories of the CBR photons, resulting 
in a reduction of the small-scale intrinsic fluctuations in CBR (e.g. Kashlinsky 1988; 
Fukushige et al. 1996 and references therein). However, this component does not cause 
additional temperature variations in CBR because of the conservation of the CBR surface 
brightness; The second one comes from the Rees-Sciama effect, which accounts for the 
changing gravitational potential during the radiation cross time, giving rise to a decrement 
of CBR temperature; And the third component corresponds to the purely relativistic 
time-delay due to the deep gravitational potential, which enables us to receive the CBR 
photons emitted at an earlier time and thus makes a positive contribution to AT/T. This 
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last component is predicted by both gravitational lensing and the Rees-Sciama effect. 

In their pioneering work, Rees & Sciama (1968) concluded that more significant 
contribution to the CBR anisotropics could be due to the time-dependent potential other 
than the time-delay effect. This has led many authors to devote their interest to the CBR 
anisotropies by the evolving gravitational potential of various nonlinear perturbations (e.g. 
Seljak 1996 and reference therein). Unfortunately, these results can hardly be applicable to 
galaxy clusters because of the unrealistic models for cluster matter distributions such as 
the the Swiss cheese model (Dyer 1976) and the "two step Vacuole" model (Rees & Sciama 
1968; Nottale 1984), while numerical simulations are limited by the dynamical resolutions 
on cluster scales especially inside core radii (Tuluie & Laguna 1995; Seljak 1996). So far, 
the only 'plausible' constraint on the CBR temperature fluctuations by the time-dependent 
potential of a rich cluster may be the work by Chodorowski (1991). He studied the effect 
using a linear potential approximation and a pure spherical infall model for cluster, and 
reached that AT/T < 6.5 x 1(T 7 . In some sense, such an estimate can be regarded as an 
upper limit on the CBR temperature fluctuations caused by the time-dependent potential 
of a cluster. 

On the other hand, numerous observations have claimed for a "settled" configuration 
of cluster matter evolution since z ~ 1. It was shown more than a decade ago that optical 
counts of clusters are consistent with no-evolution scenario for redshift out to at least 
z 0.5 (Gunn, Hoessel & Oke 1986). The same conclusion holds true for the X-ray selected 
clusters since z ~ 0.8 (e.g. Fan, Bahcall & Cen 1997; Rosati et al. 1998). Moreover, no 
significant differences in the dynamical properties have been detected between high-redshift 
and low-redshift clusters, which include the X-ray luminosity, the X-ray temperature, the 
velocity dispersion of cluster galaxies, the mass-to-light ratio, the baryon fraction, etc. 
(Carlberg et al. 1996; Mushotzky & Scharf 1997; Wu, Xue & Fang 1999 and references 
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therein). In particular, the distribution of core radii of the intracluster gas in nearby 
clusters is identical to that of distant ones (z > 0.4) (Vikhlinin et al 1998). Taking these 
observational facts as a whole, we feel that, in addition to the claim for a low- mass density 
universe, the gravitational potential of clusters is unlikely to have experienced a violent 
change since zm\. Therefore, the CBR fluctuation due to the changing or time-dependent 
gravitational potential of clusters at z < 1 could be rather small. Therefore, we need 
probably concentrate on the 'third component' and explore how large the time-delay effect 
on the CBR fluctuations by the static gravitational well of a rich cluster would be. 

The gravitational time-delay can be roughly estimated by At = 
(2GM/c 3 ) \n(AD cs D c /rl) for a pointlike mass M, where D c is the distance to M, 
D cs is the separation between M and the background source, and r is the impact 
parameter. For a rich cluster with mass of 10 15 M Q at cosmological distance and CBR as 
the background source, At ~ 10 3 - 10 4 years. Recall that the time-delay between the 
images of the gravitationally doubled quasar by a galactic lens of 10 12 M© is ~ 1 yr. A time 
delay of Ato = 10 3 - 10 4 yrs in the observer's frame corresponds to a CBR temperature 
increment of AT d /T ~ At /t ~ 1CT 7 - 1CT 6 . Indeed, such an amount of temperature 
variation is about two orders of magnitude smaller than the thermal S-Z effect. However, 
this value should be detectable with the future CBR detectors like MAP and Planck. Unlike 
the thermal S-Z effect and the kinematic S-Z effect due to the peculiar motion of galaxy 
clusters, AT d /T arising from the time-delay effect is insensitive to the impact distance, and 
thereby may be the dominant contribution to the CBR temperature fluctuations at large 
radii from the cluster centers. Alternatively, the presence of the time-delay component may 
complicate the measurement of the kinematic S-Z effect. It has been suggested that the 
thermal and kinematic S-Z effects can be separated using their different spectra especially 
at the frequency near 218 GHz where the thermal effect is zero (see Birkinshaw 1999). 
Now, one may also need to subtract properly the contribution of the time-delay effect in 
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order to extract the kinematic effect although the latter is still the dominant component. 
In particular, the time-delay effect by clusters may be comparable to the CBR perturbation 
caused by the transverse motion of the clusters as lenses (Birkinshaw & Gull 1983; Gurvits 
& Mitrofanov 1986). This will add further difficulty to the distinction between the two 
effects even if the CBR measurements can reach a sensitivity of ~ 10 _6 -10~ 7 . On the other 
hand, a quantitative analysis of the time-delay effect by clusters will be helpful for our 
estimate of various uncertainties in the CBR measurements around clusters. 



2. Time-delay effect by clusters of galaxies 

CBR anisotropy by the Rees-Sciama effect is mainly associated with nonlinear and 
strongly evolving potentials, which can be well described by the following formalism 
(Martinez- Gonzalez, Sanz & Silk 1990): 

AT 5 , , , x 2 



T 



^(0 o -0 e )-^frff-V0 + n.(^-^), (1) 



or equivalently, 



where the first term is the Sachs- Wolfe effects, the second term denotes the Rees-Sciama 
effect by a time-dependent potential well of nonstatic structure, and the third term simply 
represents the Doppler shift. We are interested in the second term generated by an isolated 
structure like a galaxy cluster. For a static potential 0(f) embedded in the expanding 
universe, we have 

f = - / : y 4 , 3 

Jv \r — r\ 

in which f and f" are connected to the comoving coordinates through r = a(t)x and 
F = a(t)x', respectively. Because of the mass conservation of Aixpif)^?, the partial 



-7- 



derivative of with respective to t reads 

9 4 = ~% (4) 
at a 

As a result, the CBR temperature anisotropy by the Rees-Sciama effect of a static potential 
4> is 

Ar Ar-^s. ( 5) 



T c 3 J a 

where the integration is performed along the light path s. If we assume that the cosmological 
term a/a remains roughly unchanged during the CBR photon cross time, i.e., the size of 
the nonlinear structure represented by the static potential is relatively small, the above 
expression can be written as 

AT a . 

At = -|y 4,ds, (7) 

where At is the relativistic time dilation due to the presence of in the framework of 
the linearized Einstein theory (e.g. Cooke & Kantowski 1975). We can also view the 
problem from a different angle: The CBR photons can be trapped in the gravitational well 
4> and separated from the expansion of the Universe for a period of At. Namely, the CBR 
photons traveling through a gravitational well can conserve their energy with respect to the 
background photons. Defining the Hubble constant as H(t) = a/a, we have from eq.(6) 

AT 

-j^ = H(t)At = VTTV c H At (8) 

where the subscript d denotes the time-delay component, and the time-delay in the 
observer's frame is Ato = (1 + z c )At with z c being the redshift of the nonlinear structure, 
Here and also hereafter we assume a flat cosmological model with Q M = 1 and H = 50 km 
s- 1 Mpc -1 . 

We now focus on the numerical computation of the time-delay effect by clusters of 
galaxies. We approximate the mass distribution of clusters by two well-known models: the 
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softened isothermal sphere model (hereafter SIS) and the cusped universal density profile 
(Navarro, Frenk & White 1995; hereafter NFW): 

a 2 1 , n . 
PSIS ~ 2^G^! ] (9) 

OcPcrit /in\ 

(r/r s )(l +r/r s ) 2 

In SIS model, a is the velocity dispersion of dark matter particles and r c is the core radius, 
while in NFW profile, p crit = 3H 2 /8irG is the critical mass density for closure, 5 C is the 
dimensionless characteristic density contrast, and r s is the scale length. To ensure the 
convergence of gravitational potential and the validity of our assumption about the limited 
effective size of nonlinear structure, we truncate the cluster at its virial radius defined by 

Air 

M(r vir ) = —r 3 mr A cPcrit , (11) 

in which A c « 200 is the overdensity of dark matter halo with respect to p cri t. The 
gravitational potentials for these two models inside follows: For SIS model 



0575 = 2<T^ 



arctamc 1 . 9 1 . 9 

+ - ln(l + x 2 ) - 1 - - ln(l + x^ r ) 



x 2 y 2 

where x = r/r c and = r vir /r c . For NFW 



(12) 



<W = -47rGr 2 p s ln(1 + :r) , (13) 

Ob 

where p s = 5 c p cr i t and x = r/r s . Inserting these derived potentials into eq.(7) and 
performing the integration along the light path across the cluster, we can obtain the CBR 
temperature fluctuations due to the static gravitational potential of clusters approximated 
by SIS model and NFW profile, respectively. 



3. Comparison 
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3.1. S-Z effect 

For an isothermal (5 model as the distribution of the hot plasma inside a cluster, the 
thermal S-Z effect is (see Rephaeli 1995; Birkinshaw 1999) 



. (l-3/3)/2 

AT T sz(#) ' " 



T, 



CBR 



= g(x v )y 



6 



(14) 



2 r 

'7* C> V / rp 



^ = ( e J_ 1)2 coth Y " 4 J ' ( 15 ) 
where i v = hu/kT C BR is the dimensionless frequency, T C br = 2.726K is the present CBR 
temperature, and 

y = 7,2 x ( I ^) ( J^) (^) , (16) 

in which n e o and T e are the central electron number density and temperature, respectively, 
and r x>c (or XjC ) is the core radius for the (3 model. The kinematic S-Z effect due to the 
peculiar motion v of cluster along the line of sight is 



= - - n 



Jcbr Vc/ 1 (-^ 



21 (1-3/3J/2 



:i7) 



There are two remarkable distinctions between the S-Z and time-delay effects: (l)Unlike 
the S-Z effect which is independent of cluster redshift, the time-delay effect varies as 
(1 + z c ) 3 / 2 ; (2)Both ATtsz and ATksz drop sharply with the outward radius from cluster 
center, while AT d is rather insensitive to cluster radius. In order to quantitatively compare 
the CBR temperature fluctuations arising from the thermal and kinematic S-Z effects and 
the time-delay effect, we take a typical rich cluster at z — 0.1, whose parameters are listed 
in Table 1, to proceed our numerical computations. For the NFW profile, we adopt the 
typical values of a(= 4irGfj,m p p s r 2 /kT e ) and r s found from the fitting of the NFW expected 
X-ray surface brightness profiles of clusters to the observed ones (e.g. Ettori & Fabian 1999; 
Wu & Xue 2000). The resultant CBR temperature variations AT(#)/Tcbr are shown in 
Fig.l for two different choices of the truncated cluster radii: r = r vir and r = 10r„j r . 
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EDITOR: PLACE FIGURE |T| HERE. 

It appears that for a typical cluster and an observing frequency v = 32 GHz, the 
orders of magnitude of the maximum CBR temperature fluctuations by the thermal and 
kinematic S-Z effects and the time-delay effect are, respectively, 10~ 4 , 10~ 5 and 10 -6 . 
Although a temperature fluctuation of as low as 10 -6 will be detectable with the future 
space experiments like MAP and PLANCK, the fact that AT^/Tcbr is a slowly varying 
function of radius puts the actual measurement of the time-delay effect into a difficult 
position. One possible way is to measure the temperature difference between two points 
separated by an angle A9. The signature of \ATd(0 — A9/2) — ATd{6 + A9/2)\ could be 
identified if the S-Z effects can be removed. Yet, the kinematic S-Z and time-delay effects 
may become to be indistinguishable as a result of the frequency-independent property, 
unless the CBR temperature profiles with a sensitivity of at least 10~ 7 can be obtained. In 
other words, the presence of the time-delay effect may yield an uncertainty of ~ 10% in the 
measurement of the central kinematic S-Z effect of clusters. 



Table 1. Cluster parameters 



a (km/s) 


P 


r c (Mpc) 


n e0 (cm 3 ) 


T e (keV) 


v (km/s) 


a r s (Mpc) 


1100 


2/3 


0.25 


3 x 10~ 3 


7 


500 


10 0.8 
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3.2. Transverse motion of clusters 

In the Rayleigh- Jeans limit, the magnitude of the temperature fluctuation AT v /T C br 
due to the transverse motion of a cluster with velocity v can be estimated through 
(Birkinshaw k Gull 1983; Gurvits k Mitrofanov 1986) 

where 5(9) is the deflection angle produced by the projected cluster mass within 9 along the 
line of sight. For the SIS and NFW models, we have 



fa 2 \ J9 2 + 9 2 C -9 C 
Ssis(e) = 4n i-U f , (19) 



and 

„2 



we) = =. mx , «• SF* • ' (20) 



respectively, where S is the angle distance of r s . AT 1 ,(6 i )/Tcbr produces a two-sided pattern 



around the moving cluster. Unlike the time-delay and S-Z effects, CBR will be unaffected 
by the transverse motion of the cluster if we look through its central region because of 
5(0) = 0. The maximum amplitude of AT V /T CBR relative to AT„(0)/T CB r for a typical rich 
cluster is 

|| KO ' 9xlo "(lo^)(lo%^) 2 - SIS; (21 > 

&T — » U x 10- f^^) NFW. (22) 



T CB R Vl0 3 kms-V VlO/ VTkeV, 

Apart from their very different CBR patterns, AT KS z/T C br, AT v /T C br and AT d /T CBK 
are all achromatic, and the latter two effects also have the same order of magnitude. 
Consequently, the uncertainty in the measurement of ATksz due to the combined effect of 
the time-delay and the transverse motion of clusters can become even larger than ~ 10%. It 
was suggested (Birkinshaw k Gull 1983; Gurvits k Mitrofanov 1986) that the detection of 



AT v /T C br can be used as a method for measuring the peculiar velocities of clusters. Such 
a motivation can now be complicated by the time-delay effect unless the detailed patterns 
of the CBR anisotropics around clusters can be well mapped. 



4. Discussion and conclusions 

Indeed, the CBR temperature fluctuation caused by static gravitational potential of 
a rich cluster is very small, AT d (6)/T CBR ~ 10 -6 -10 _7 , which is of 2-3 (1-2) orders of 
magnitude lower than the thermal (kinematic) S-Z effect, but is nevertheless comparable 
to the effect produced by a transversely moving cluster as the gravitational lens. The 
signals of the time-delay effect and the transverse motion of clusters may remain to be 
indistinguishable from the kinematic S-Z effect in current S-Z measurement, unless one 
can acquire the detailed CBR temperature profile across clusters with a sufficiently high 
sensitivity of ~ 1CT 7 . The presence of the time-delay effect and the transverse motion of 
clusters may lead to an uncertainty of ~ 10% in the measurement of the kinematic S-Z 
effect due to the peculiar motion of clusters along the line of sight, which gives a sense 
of how accurate and robust one can use the kinematic S-Z measurement of clusters for 
cosmological purpose. 

So far, we have not included the contribution from the time-dependent potential of 
clusters. This is mainly based on the recent observations that the significant evolution of 
dynamical properties of clusters has not been found since z ~ 1. Therefore, our conclusion 
may not hold exactly true if clusters are still in the process of formation where the free infall 
plays a dominant role. Our derived CBR temperature fluctuation due to the time-delay 
effect can be comparable to that produced by the time-dependent potential of clusters (e.g. 
Chodorowski 1991; Tuluie & Laguna 1995; Tuluie, Laguna & Anninos 1996). Nevertheless, 
in a similar way to the S-Z effect, the AT/Tqbr caused by the changing gravitational 
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potential of a cluster shows a sharp drop along outward radius. So, it would be possible to 
isolate the time-delay effect from the measured CBR fluctuations behind a cluster if one 
can have the high-sensitivity CBR temperature profile. 

Because of the unique property of the 'long-distance' effect, i.e., AT d /T C BR depends 
approximately on ln(l/0), one may worry about the issue of whether the time-delay effect 
by clusters can contaminate the global CBR power spectrum measured at small angles 
(~ 10-100 arcminutes) because clusters are a rare population in the Universe. Without 
a sophisticated statistical study of the effect on the CBR spectrum, our computation in 
the present paper suggests that the time-delay effect by clusters is unlikely to produce a 
noticeable contribution of as high as few times 10~ 6 to the CBR anisotropics, which is 
compatible with the Rees-Sciama effect generated by large-scale matter inhomogeneities 
according to numerical simulations (Tuluie & Laguna 1995; Tuluie et al. 1996; Seljak 1996). 
Nevertheless, future precise CBR temperature measurements on smaller angular scales 
should allow the time-delay effect of clusters to be included. 

We gratefully acknowledge Tzihong Chiueh for useful discussion, and the science editor 
of the journal, Prof. R. T. Vishniac, and an anonymous referee for valuable comments and 
suggestions. This work was supported by the National Science Foundation of China, under 
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Fig. 1. — Radial CBR temperature variations generated by the thermal (dotted lines) and 
kinematic (dashed lines) S-Z effects and the time-delay (solid line) effect for a cluster at 
z = 0.1. For the thermal S-Z effect an observing frequency of v = 32 GHz is assumed. 
Cluster properties are summarized in Table 1. Dependence of the effects on the truncated 
radii is shown for r cut = r vir (upper panels) and r cut = I0r vir (lower panels), respectively. 



